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Abstract. - Charge transfer statistics of quantum particles is obtained by analysing the time 
evolution of the many-body wave function. Exploiting properly chosen gauge transformations, 
we construct the probabilities for transfers of a discrete number of particles. Generally, the 
derived formula for counting statistics differs from the one previously obtained by Levitov 
et al. (J. of Math. Phys. 37, 4845 (1996)). The two formulae agree only if the initial 
state is prohibited from being a superposition of different charge states. Their difference is 
illustrated for cases of a single particle and a tunnel junction, and the role of charge coherence 
is demonstrated. 



Recently the question of counting statistics of charge transfer has attracted considerable 
interest due to its relevance to electronic transport in nanostructures where the discreteness of 
the electronic charge is reflected in quantum transport properties 'l"2"3"4'5'fc)' . Inspired by the 
concept of photon counting in optics, counting statistics of particles addresses a fundamental 
question of quantum transport, viz. the probability distribution for the number of charges 
transferred between different spatial regions of a system in a given time span. The objective 
is to get complete information about the fluctuations in particle currents, i.e., correlations of 
any order. Generally, cross-correlations in charge transfers in different conducting channels of 
a mesoscopic system are of interest (Tj. In a seminal paper by Levitov, Lee and Lesovik PP, a 
formula for counting statistics was proposed by considering a gedanken experiment in which a 
spin is coupled to the electrons in a quantum wire whose transfer statistics are to be counted. 
The precession of the spin then counts the number of electrons passing either to the left or 
to the right of a chosen point in the wire. Applying quantum measurement arguments to a 
system interacting with an idealized measuring device, counting statistics was also considered 
by Nazarov and Kinderniann 'W . 

Here, we shall develop counting statistics from a different point of view; instead of analysing 
a measurement process, we extract information about particle transfer directly from the wave 
function of a many-body system. We derive a formula for counting statistics which turns out 
to differ from the one of Ref. jlj . Our approach allows us to establish the circumstances under 
which the counting formula obtained in Ref. P is not applicable. 

In classical mechanics, the notion of counting statistics is unproblematic: when complete 
information about a system is known - the trajectories of individual particles are known - 
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there is a unique answer to the question of how many particles are transferred from, say, 
the left to the right in a given time interval. In quantum mechanics the situation is not 
so innocent. Even when full information about a many-body system is available, i.e., its 
time-dependent wave function is known, there is no straightforward algorithm to extract the 
probabilities in question, since there is no quantum operator representing the number of 
transferred particles. To circumvent this difhculty, a gedanken experiment was in Ref. 1 used 
as the basis for establishing a counting formula: In the experiment, the rotation x(A) of a spin 
coupled to the charge current via a gauge field was "measured" as a function of the coupling 
constant A (and measuring time interval t). Assuming the generating function x{^) to be 
27r-periodic, the Fourier coefficients in front of exp(iTOA) were interpreted as the probabilities 
for the passage of m particles. We point out that the interpretation of the experiment is based 
on ingenious intuition rather than following unequivocally from the principles of quantum 
mechanics. Besides, it contains an ambiguity: as we show later, x{^) may have exp(±iA/2) 
components, which could then be interpreted as half-integer charge transfers. Also, positive 
definiteness of the prescribed probabilities cannot be established. A gedanken experiment, not 
being a realistic one, is in fact a vehicle for analysing the wave function of a system. Therefore 
we shall develop an approach to counting statistics of quantum particles which is based solely 
on an analysis of the wave function. 

We consider a system partitioned into two parts by the plane at x = 0, referred to as left 
and right. We want to tag the particles with a "non-demolishing marker," i.e., a marker that 
does not disturb the quantum dynamics. The marker should provide information on whether 
a particle has crossed the interface and in what direction. As explained below, this can be 
achieved by introducing the gauge transformation 



U\ = exp 



-Xk] 



(1) 



where X]^ is the coordinate of the fc-th particle of the many-body system. In order to demon- 
strate how the gauge transformation serves as a marker, we consider the case of a single 
particle subject to a potential. Let 4'l{x) and iPr{x) denote normalized initial wave packet 
states located only on the left or right, respectively. A state is evolved in time according 
to ?/'(t) — £V(0), where £ — exp[— iTir] and TL is the Hamiltonian of the system. For each 
initial state, specified at time t = Q, the time evolution operator £ produces a state which is 
a coherent superposition of left and right components: 

£'>pL = ^L^L + fpL^R , £^R = ll^R~.L + fpR^R , 

where the last symbol in the subscript on the r.h.s. indicates the location of the wave packet 
to be on the left (L) or right (R). Of importance will be the gauge transformed evolution 
operator 

£x = ul£Ux. (2) 
Letting it operate on the considered initial states, the following marked final states emerges 

£\ll^L = V'L^L + e*-^ i^L^R , £xll^R = e~*^ i^R^L + iJ-'R^R , 

indeed states exhibiting the intended transfer marking. One immediately realises that the 
weights IIV'l-.lIP = and Wi/jl-^rW^ = {i^L~>R.\ipL~*R) are the probabilities for 

the charge transfers m — and m — 1, respectively. Analogously, ||V'_r^l||^ and HV'fl^ijlP 
are the probabilities for transfers m — —1 and m = for the initial state on the right. One 
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is able to extract the components, V'm , of the final states E^^l and £tpR which corresponds to 
m particle transfers from left to right by the following operation 

= /ge-™^fA^W, (3) 



where V'''^-' is the initial state, i.e., either ipL or ijjR. The probabilities for the possible particle 
transfers can therefore also be expressed on the form Pm = (V'mlV'm)- 

Next we demonstrate that the procedure works for an arbitrary initial state and consider 
a superposition of the wave functions located on the left and right, V^^^^ = AipL + The 
operation in Eq. Q then produces the three charge transfer states 

AipL^ii , m = 1 ; 
V'm = { Ai}L^L + Bipii^ii , m = ; 

BipR^L , m=-l. 

We note that the weight of these states, H'f/'mlP — (V'mlV'm), are the probabilities for transfers 
of TO particles for a general initial state. Indeed, = is the product 

of the probability initially to be on the left side and the conditional probability to transfer 
from the left to the right, and analogously for the terms (V'olV'o) and 

Using Eq. ©, the transfer probability, P„i (V'mlV'in); can be expressed as 

P„. = /^|e--(^(")|4_.£,,.|^(0)>, (4) 



and the generating function, XtW = ^ PmS^"^^ , reads 

m 

Xr(A) = 1 ^ Xr(A, A) , Xr(A, A) = {£l_^£A+i) , (5) 



where the averaging is with respect to the initial state -;/''•*'•' , or the density matrix of the system. 
This expression is valid for a many-body system provided the corresponding gauge-transformed 
evolution operator Eq. |2Jl is used. Inverting the argument, Fourier transformation with 
respect to A of the 27r— periodic function Xr(A), Eq. (Q, generates the coefficients Pm (for 
integer to's) which are guaranteed to be positive, and ^™ ^ ^- Their meaning is that of 
probabilities for integer charge transfers. 

Expressing the evolution operator via the Hamiltonian, the integrand of the generating 
fimction in Eq. ((SJ can be written as 



Xr(A, A) = ( Ta'Cxp 



i / dt' n-,(^t'){t') 



(6) 



where Ti.^ — U^TiU^, and H is the Hamiltonian for the system, and the Keldysh contour Cr 
proceeds from t- = —oo to r and back again as from r to — oo; Tr- denotes the time 
ordering on the contour. The projecting gauge field 7 is zero outside the measuring interval, 
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and 7(t^) = A ± -I for < t' < r, and the average in Eq. (0 is taken with respect to the 
density matrix in the far past. The expression in Eq. © can be evahiated using standard 
field theoretical methods. 

The formula, Eq. |^ , differs from the generating function proposed by Levitov et al. [T]|^ . 
The latter, denoting it x^, is obtained by setting A = in Eq. © 

x\{\) = {£\S.). (7) 

2 2 

To verify that the two counting formulas are not equivalent, we calculate 'X^{\) for the single 
particle case where the initial state is the previously considered superposition of right and left 
located wave packets, and obtain 

X^^(A) = xr (A) + 4z sin (A*B {^L^R\i^R^R)) ■ (8) 

Indeed, Xr differs from our generating function. We find that the difference is the additional 
term in Xr which is 47r— periodic in A. In Ref. 1 where the procedure for charge transfer 
counting were based on the Fourier expansion of ^{X), this implies that half- integer charge 
transfers would occur. 

To investigate further the difference between the two approaches, we consider counting 
statistics from a different perspective. Introducing the Hermitian operators 

Pn= f^e-'^^U^ , n = 0,±l,..., 



through the marker gauge transformation, Eq. we realize their meaning by noting that 
the operator V,i projects a state {ip) onto the component lipn) = 'Pnl'ip) which corresponds to 
exactly n particles on the left. These projection operators, similar to the ones introduced by 
P. W. Anderson in superconductivity, turn out to be suitable tools for the kind of vivisection 
of a quantum state needed to obtain the probability distribution for discrete charge transfers. 
Being states with definite particle number on the left, the projections IV'n) are eigenf unctions 
of the operator U^, C/^IV'™) = e*'''"'|^„). This property can be expressed on operator form as 
Vn = e~'^^'^UjVn, whereby J2n e^^Vn = Uj . Consequently, Eq. can be transformed into 

Prn = 5](V'"'|n.f^K-™^^K|V'(°^), (9) 
n 

producing a different way of expressing the transfer probability. According to quantum me- 
chanics, the matrix element {'ip^°^Vn£'''Pn-m£'Pn\ip^°'') = \ \'Pn-m£'Pntp^°'> is the probability 
for the transition from a state with n particles on the left to a state with n — m particles on 
the left. The quantity Pm is thus the probability for a transfer of m particles to the right 
in a time span r given that a measurement of the charge state is performed initially. For an 
arbitrary mixture of states, 

£Vn) , (10) 

n 

where, as in Eq. (|SJ, the average means taking trace with respect to the density matrix po at 
time t = when the counting is initiated. (We recall that the evolution operator £ evolves 
the system from time t — Q to t = t.) For a classical statistical ensemble, this is how the 
statistics of particle transfers is evaluated, and we conclude that Eq. (|10|) and therefore the 
generating function, Eq. lO, indeed has the correct classical limit. 
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In terms of the charge projection operators, the generating function Eq. ((SJ becomes 

xr(A) = V (n.^l.^Ano- (11) 

n 

We infer from this expression that if po is diagonal in the representation of the charge states, 
one of the projection operators can be removed from Eq. and the sum over the remaining 
projectors is unity; equivalently, Xt{^t^) does not depend on A and the integration with 
respect to A can be omitted. In this case, the generating function, Eq. (jlll) . reduces to the 
form on the r.h.s. of Eq. {T)), which is identical to the result of Ref. pp. The physical origin 
of the difference between the two formulae is thus the charge off-diagonal components of the 
density matrix. As one can show, it is the latter which produces the unphysical 47r-pcriodic 
part of the generating function of Ref. P . 

We observe that the simple physical picture in which the current is build of transfers 
of integer number of particles meets with difficulties for a general initial state, when the 
system may be in a superposition of different charge states, and po is non-diagonal in charge 
space. Indeed, one expects that Q{t) — e rnP^n should equal the average charge transfer, 

{I{t))dt, I being the Heisenberg current operator. Calculating Q from Eq. ^ or Eq. 
Q{t) — '}2n{^riI{t)Vn) ■ Clearly, Q is identical to (/) only when po is charge diagonal. A 
similar difhculty emerges for the generating function of Ref. ^j: the expression Q — '^^m "^^m 
gives the correct expectation value for the transferred charge only if one allows m to assume 
half-integer as well as integer values, the former being due to the 47r-periodic part of Xr (-^) 
generated by the charge off-diagonal elements of po- The analysis points to an ambiguity in 
counting statistics, a trade-off between having a probability distribution for discrete charge 
transfers and the generation of proper current correlation functions. The resolution is shown 
to be tied to the charge structure of the initial state of the system. 

As an illustration, we evaluate the counting statistics for a tunnel junction using Eq. (O; 
the problem was considered previously in Ref. 9 using the generating function in Eq. Q . The 
system consists of two weakly connected metallic regions, left and right. The Hamiltonian 
reads H = Ho + V'^ , where Ho refers to the isolated regions and V'^ = Vr^i + Vi^r is 
the tunneling part, where Vr^i ( V/wr) describe transitions from left to right (right to left), 
Vr^i — V'j|_j,. The gauge transformation Uj affects only the tunneling part, and in Eq. lO, 

n^^Ho + , = e'^Vr^i + e-'-'Vi^r ■ (12) 

As in Ref. [H], we consider only the leading contributions with respect to tunneling in Eq. (jSJ, 
and W{X,A) = lnx(A, A) can be evaluated as 

W{X,A) - /tk J diidt2l>^t,)(ti)l>^t,)(t2)\ , (13) 

^ — oo ' 

where V'^{t) is the tunneling operator in the interaction picture; for given A and A, the 
projecting field 7 is the function of the Keldysh time t± introduced below Eq. @. One 
obtains two contributions, W = Wmm + Wmp- The first term, Wmm, originates from the 
time domain where both time arguments ti and t2 in Eq. (|13|l belong to the measuring interval 
from to T, and the second one, Wmp, is the contribution from the region where one of the 
times ti or t2 is prior to the start of the measurement, t = 0. 

The term Wmm is 27r-periodic in A and does not depend on A. It coincides with the result 
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of Ref. 0: Wmm = (e'^ - 1) w+{t) + (e"*^ - 1) w_(r) where 

e'" - i ^ 

2ttt 



w+{t) = 2^T y" dEdeTE_,E+ n'E_ (l - n^^) A,(e) , A,(e) - 



Here E± — E zt ^, and are the electron distribution functions for the energy E in the 
left and right regions respectively, and 

te,e' = J2 - - v), 

p,p' 

£p being the single particle energy; one obtains W- by substituting n''^ for (1 — n}'^) in w+. 
The second term, Wmp, has the form 

WMP = 2ism^^{e'^WMp{T)) , 



where 



dEdeTE_^E+ [nE_ -riE^j ^^^^^ , e+ = e + iO 



and Wmp is thus a 47r— periodic function of A. 

The generating function factorizes, x — Xmm Xmp, where xmm = e^'^'" is identical to 
the result in Ref. [H], and 



Xmp 



(r,A)= ^'"'^Jli^MPir)) , (14) 



Jm being the Bessel function. If Te^e' is featureless on the scale of the Fermi energy, Ep, 
Wmp is a constant once r ^ h/Ep, 

WMp{oo) = , (15) 

Jrij' 

where Rt is the tunneling resistance and Rq = 27r?i/e^. 

The calculation for the tunnel junction qualitatively agrees with the result we obtained in 
the single-particle case. We observe again that the formula of Ref. UEIj which is identical to 
Xr(A,A = 0), contains 47r-periodic terms. According to the derivation, they originate from 
tunneling events which occur before the measurement started and create charge off-diagonal 
elements in the density matrix by the time t = 0. For large enough measuring times, Wmp 
saturates unlike Wmai which grows linearly with time, and Wmp represents memory of the 
initial state of the system, amounting to intrinsic voltage independent charge fluctuations. 
The latter, however, need not be small, as seen from Eq. H15|l . According to our analysis, the 
results in [H], based on Ref. pj, are valid only if the two electrodes are not connected before 
the measurement whereby charge superposition is prevented. 

In this paper we have reconsidered counting statistics applying the rules of quantum me- 
chanics. We confirm the counting formula of Levitov et ai, but only for cases where the initial 
state of the system is charge diagonal, i.e., when superposition of different charge states is ab- 
sent. Our approach leads to a novel formula for the probability distribution of integer charge 
transfers which is valid for a general charge coherent state. 

The role of charge coherence, which is a ubiquitous feature of any many-body quantum 
state, should be examined in each particular case. For a system of non-interacting particles. 
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one may argue that details of the initial state are of minor importance since the net contri- 
bution of the non-diagonal elements tends to average out. Our tunnel junction results show 
how it comes about in this particular case: Even though charge coherence is present in the 
initial state at t = 0, created by prior tunneling events, its contribution expressed by wmp 
does not grow with time, diminishing in importance at large measuring times. Nevertheless, 
the charge coherence present at the start of counting, does noticeably change the statistics 
in accordance with our formula Eq. (|14|l . especially for short measuring times. For interact- 
ing systems the situation is similar, provided the charge structure of the quantum state can 
be expressed in terms of quasi-particles. An important counterexample is a superconductor, 
where the superposition of different charge states is rigidly maintained, as required by the 
number-phase uncertainty relation. Although the theoretical objects entering the formula 
of Levitov et al. can be calculated for superconductors 10 and even measured jDEli they 
cannot be interpreted as charge transfer probabilities. Since charge off-diagonal elements are 
important, the formula of Ref. I, cannot be used to calculate the statistics of charge transfer 
related to the current of Cooper pairs. The fact that the previous counting formula leads to 
negative "probabilities" in the case of a Josephson junction is understandable in view of 
this observation. 

In conclusion, we have shown that by using gauge and charge projection operators to 
analyze the structure of a quantum state of an arbitrary system, one is able to construct 
a probability distribution for charge transfers of particles obeying quantum dynamics. The 
constructed function is a proper probability distribution, i.e., positive definite and normalized, 
and the probability distribution for charge transfer counting of classical mechanics emerges 
in the correspondence limit. The charge transfer in a tunnel junction is considered, and the 
modification of counting statistics due to charge coherence has been demonstrated. 
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